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f^ ■ We analyse the phase diagram of a quantum mean spherical model 

£h ■ in terms of the temperature T, a quantum parameter g, and the ratio 

p = —J2/J1, where J\ > refers to ferromagnetic interactions between 
first-neighbour sites along the d directions of a hypercubic lattice, and 
J2 < is associated with competing antiferromagnetic interactions 
between second neighbours along m < d directions. We regain a num- 
C"} ' ber of known results for the classical version of this model, includine; 

£— n ■ the topology of the critical line in the g = space, with a Lifshitz 

■<sj- , point at p = 1/4, for d > 2, and closed-form expressions for the decay 

(T) I of the pair correlations in one dimension. In the T = phase diagram, 

there is a critical border, g c = g c (p) for d > 2, with a singularity at 
the Lifshitz point if d < (m + 4)/2. We also establish upper and lower 
critical dimensions, and analyse the quantum critical behavior in the 
neighborhood of p = 1/4. 



1 Introduction 

The spherical model of magnetism has been used as an excellent laboratory to 
test ideas and concepts of phase transitions and critical phenomena [lj[2j[3j. 
There are several versions of the original model, including proposals of a 
quantum spherical model to correct some of the unphysical results at low 



temperatures [I] [2| [5] [I] • The effects of frustration [S][H], random fields [TO] . 
and of disordered exchange interactions [H][I2], have also been analysed 
in the context of quantum spherical models. With a view to describe the 
crossover between classical and quantum critical behaviour, Vojta [6] used a 
standard scheme of canonical quantization to analyse a quantum version of 
the ferromagnetic mean spherical model. We were then motivated to revisit 
this problem, with the addition of competing ferro and antiferromagnetic 
interactions, and the perspective to analyse a quantum Lifshitz point. 

The mean spherical model, which has been originally proposed by Lewis 
and Wannier [2j, is given by the partition function 
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where /3 = 1/ (ksT), T is the temperature and ks is the Boltzmann constant, 
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/i is a suitable chemical potential, I is a lattice vector, and {S-f} is a set 
of continuous spin variables running over the N d sites of a <i-dimensional 
hypercubic lattice. The model Hamiltonian is written as 
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where ( k , I ) labels a pair of lattice sites, the exchange parameter J 
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t-f 



k, I 



depends on the distance between sites k and / , and H is an 



external field. In this formulation, the chemical potential /i comes from the 
mean spherical condition, 
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and it is well known that exact solutions for the thermodynamic functions 
can be obtained from the standard diagonalization of a quadratic form [3] . 
In a quantum version of this mean spherical model [1] [B] , the spin variable 

S-f becomes a position operator at lattice site / , canonically conjugate to a 
momentum operator P-f, with the commutation relations 
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where 8-? -7+ is a Kronecker delta and we assume that H = 1. We then add a 

< , k 

term of kinetic energy, depending on a quantum parameter g, and write the 
quantum quadratic form 

^=^E P ?- E J t,?% 5 T-#E 5 ?+^E4' ( 5 ) 
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which can be diagonalized by a canonical method [5J, leading to a solution 
of the problem for a general ferromagnetic pair interaction. At finite tem- 
peratures, the critical behaviour is essentially unchanged with respect to the 
classical spherical model. At zero temperature, depending on the parame- 
ter g, there is a quantum phase transition characterized by new (quantum) 
critical exponents. Also, the introduction of quantum fluctuations leads to a 
correction of the unphysical behaviour of the entropy at low temperatures. 

We report an analysis of this version of the quantum mean spherical 
model in the presence of competing interactions. We consider ferromagnetic 
interactions, J± > 0, between pairs of first-neighbour sites along the d di- 
rections of a hypercubic lattice, and antiferromagnetic interactions, Ji < 0, 
between second-neighbour sites along m < d directions. Classical versions of 
this model [13] [H] [15] [16] , as well as more elaborate mean spherical models 
with competing interactions [17] . have been studied by several authors. For 
m — 1, we regain a spherical analogue of the Axial- Next-Nearest-neighbour 
Ising, or ANNNI, model [IH][IH], which is known to display a rich phase 
diagram, including a Lifshitz point, in terms of the temperature T and a 
parameter p = —J2/J1 that gauges the strength of the competing interac- 
tions. We then analyse the T — p — g phase diagram, for different values of 
m, in particular the T = behaviour, and establish the critical dimensions 
and critical exponents associated with this quantum model system. In the 
classical case, g = 0, we confirm a number of results, including a singularity 
of the critical border at the Lifshitz point for 2 < d < (m + 6) /2. In one 
dimension, we derive analytic expressions for the decay of pair correlations, 
and determine the region of modulated behaviour in the T—p phase diagram. 



2 The quantum mean spherical model with 
competing interactions 

This problem can be treated either by a conventional reduction to a system 
of coupled harmonic oscillators or by a judicious application of the method 
of path integrals [5] [20] [21]. Let us first use the representation in terms of 
harmonic oscillators. We then introduce bosonic operators a] and a/ to write 
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where we have omitted the vector notation. We now assume periodic bound- 
ary conditions, and change to a Fourier representation, 
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where a„ and a\ are bosonic operators, 
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and the sum is over the d-dimensional vectors q belonging to the first Bril- 
louin zone. We then write the quantum quadratic form (|3j) in the Fourier 
representation, 
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where the sum runs over all lattice vectors. The final diagonalization of this 
quadratic form comes from the introduction of new bosonic operators, c q and 
cL according to a well-known Bogoliubov transformation. We then have 
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which requires that the chemical potential // should be larger than a certain 
limiting critical value /i c , 
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In analogy with a system of harmonic oscillators, we then write the par- 
tition function 



Z n (P,H,/j) = exp 
and the free energy per site, 
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with /i > yU c = J(q c )/2. The spherical constraint, given by Eq. ([3]), from 
which we determine the chemical potential /i, is written as 
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In the classical limit, g — > 0, we regain most of the well-known results for the 
mean spherical model. The classical limit of this quantum free energy, how- 
ever, includes an extra term of the form ln(/3g), which corrects the classical 
behaviour at low temperatures. The particular limit J(q) — > corresponds 



to free quantum (spherical) rotors, with a finite energy gap, in contrast to 
the usual Heisenberg-Dirac spins. 

These expressions also come from a straightforward application of the 
path integral formalism, which has been widely used to treat quantum sta- 
tistical problems [5J[20J. It is then interesting to write the Lagrangian asso- 
ciated with this problem, 

£ = ^E^ + ^E^a^+^E^-^E^ 2 - ( 18 ) 
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In the imaginary time formalism, with t — > —ir, the partition function is 
written as 
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where the first integral includes periodic conditions, 5/(0) = S[(f3). We now 
introduce the Fourier transformation, 



Si(r) 



^ exp [i(ql + tw)] S q 



w 



(20) 



q,w 



where the vector q belongs to the first Brillouin zone, and w = w n = (2wr)//3, 
with integer n, is a Matsubara frequency. We then write the partition func- 
tion 
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where J(q) is given by Eq. (TTT1) . and we have omitted the vector notation. 
If we calculate the Gaussian integrals, and use the identity 
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it is straightforward to regain the partition function given by Eq. (1151) . 

This general solution works for all forms of distance- dependent inter- 
actions, J-g -f = J ( k — I ) . We now consider ferromagnetic interac- 
tions, Ji > 0, between pairs of first-neighbour sites along the d directions 
of a hypercubic lattice, and antiferromagnetic interactions, J<i < 0, between 
second-neighbour sites along m < d directions. The Fourier transform of the 
exchange interactions is given by 
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J{ q ) = 2 J\ 2_. cos Qj + 2 J2 2_, cos 2<?i 
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where q = (gi, g 2 , ..., g^) is a wave vector in suitable (dimensionless) units. 
The maximum of J(q) depends on the ratio p = —JijJx- If P < 1/4, the 
maximum is located at the critical value q c = 0, as in the simple ferro- 
magnetic case. If p > 1/4, the maximum of J( q ) is given by the critical 
vector 

~tc = (q c i, Qc2, -,q cm , 0, ...0) , (24) 

where 
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The special case p = 1/4 corresponds to a Lifshitz point of degree m. 



3 Phase diagrams and critical behaviour 

In zero field, H = 0, the paramagnetic critical boundary in the T — p — g 
space comes from the spherical constraint, given by Eq. ffT7|) . supplemented 
by the critical limit of the chemical potential, fi = /i c , given by Eq. ( FT4"|) . We 
then write 
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where 



w c (q) = g 1/2 [J(f c ) - J(f ) 
In the classical limit, g — > 0, we have 
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from which we obtain the critical temperature as function of p, for all values 

of d and m, 
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for p < 1/4, and 
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for p > 1/4. 

There is a long history associated with the calculations of similar lattice 
Green functions [3] [22j . From the identity 
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where n is an integer, and a > 0, we write integral representations for 
I(p,d,m), which are convenient to carry out an asymptotic analysis. For 
d < 2, the divergence of these integrals indicate that T c = for all values of 
77i and p ytz 1/4. Also, we have T c > for d > 2, and for all values of m and 
p ^ 1/4. In particular, at the Lifshitz point, T c > for d > (m + 4)/2. We 
now consider the graphs of T c = T c (p) versus the parameter p. It is easy to 
write an expression for dT c /dp, for p < 1/4 and p > 1/4, and to show that 
there is common tangent at the Lifshitz point, p = 1/4. For example, for 
d = 3 and m = 1, we have 
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Figure 1: Classical T — p phase diagrams, for dimensions d — 3 and d — 4, 
near a Lifshitz point (p = 1/4), with m — 1. The critical line separates 
ordered and disordered phases. Along the critical line, q c = for p < 1/4, and 
q c 7^ for p > 1/4. The inset shows a magnification in order to emphasize 
the singular behavior of the paramagnetic border near the Lifshitz point in 
three-dimensions. 



We then use the identity fl32l) . with n — 2, to write an integral representation 
from which it is easy to show that the common tangent at the Lifshitz point, 
p = 1/4, is infinite for (i) d < 2, and (ii) d > 2, with d < (m + 6)/2, 
which includes the analogue of the ANNNI model (d = 3 and m = 1). In 
the numerically obtained graphs of figure 1, we sketch typical profiles of the 
critical line T c = T c (p) for dimensions d = 3 and d = 4, and m = 1. Note 
the smooth behaviour of this paramagnetic border for d — 4 (and m = 1). 
The scale of this figure, however, is not enough to show the sharp singularity 
at the Lifshitz point for d = 3 (and m — 1), as pointed out in a sketch by 
Hornreich [23] . 

In the zero-temperature limit, T — > 0, Eq. ( )26l) can be written as 
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from which we obtain the critical quantum parameter, g c , as a function of p, 
for all values of d and m, 
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for p > 1/4. We now use an analytic continuation of the identity ( l32l) . for 
non integer values of n. From a similar analysis of convergence of these 
expressions, it is easy to show that there is a quantum phase transition 
(g c t^ 0) for d > 2, independent of the value of m. In particular, there is a 
common derivative dg c /dp at the Lifshitz point, p = 1/4, with a singularity 
for 2 < d < (m + 4)/2 (and a smooth behaviour for d = 3 and m — 1). In 
figure 2, we sketch typical profiles of this critical line in the g — p plane for 
dimensions d = 2 and d = 3, and m = 1. 

The critical behaviour in zero field, H = 0, comes from an asymptotic 
analysis of the spherical constraint in the neighbourhood of the transition. 
At finite temperatures, T ^ 0, in the limit \x — > J(q)/2, we have 
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At zero temperature, T = 0, we have 
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We now expand J(q) as a Taylor series about q = q c , in the classical and 
quantum cases. Although q c depends on the parameter p, the convergence 
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Figure 2: Quantum phase diagram (T = 0) near a Lifshitz point with di- 
mensions d = 2 and d — 3, and m — 1. Along the critical line, <7 C = g c {p), 
we have g c = for p < 1/4, and ^ 7^ for p > 1/4. The inset shows a 
magnification of the border near p = 1/4 in two dimensions. 
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of the sums in the right-hand side of Eqs. (157)) and (155)) does not depend 
on p, for p 7^ 1/4. For finite temperatures, T ^ 0, the sum converges if 
d > 2 (which determines the lower critical dimension of the classical case), 
regardless of the value of m. At T = 0, the sum converges for d > 1, which 
leads to the lower critical dimension of the quantum case. 
Let us consider some special situations. 

3.1 Critical behaviour at finite temperatures and p ^ 

1/4 

In zero field, H = 0, at finite temperatures, T^O, and for p ^ 1/4, we use 
Eq. (1571) and perform a Taylor expansion about T = T c and p = p c . We 
then have the asymptotic expression 



« ^ ^ / d d g :— -U- p c ) T c I d d q ? . (39) 

We now expand the integrands about q = q c , and obtain 

{r, for d > 4 
&, forrf = 4 , (40) 

r 2 , for d = 3 

where t = (T — T c )/T c gives the distance from the classical critical point. 
From these equations, supplemented by standard scaling considerations, it is 
possible to calculate the usual critical exponents associated with the classical 
spherical model. Although q c assumes different values for p < 1/4 and 
p > 1/4, the asymptotic behaviour is the same, regardless of the value of m. 

3.2 Critical behaviour at T = and p ^ 1/4 

In analogy with the calculations for finite temperatures, we write 

o*(g-9c)f d d q- 1 -— J - 2 -{n- » c )g c f d d q- l ——. (41) 

J {f^-¥) J {^-¥) 

From an expansion about q = q c , we obtain the asymptotic behaviour 

5, for d > 3 
fa-0c)~< £s, ford = 3 , (42) 

5 2 , for d = 2 
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where 5 = (g — g c )/g c gives the distance from the quantum critical point. As 
in the case of finite temperatures, Eq. ( |42p holds for p ^ 1/4 and any value of 
m < d. It is easy to use scaling arguments in order to obtain the (quantum) 
critical exponents for p ^ 1/4. With the necessary reinterpretations, and 
although critical dimensions are different, these values are in agreement with 
results of Vojta for the quantum ferromagnetic case [6]. 



3.3 Critical behaviour for p = 1/4 

At the Lifshitz point, p = 1/4, the maximum of J (q) is still given by 
q c = (0, 0, 0, ■ • • ,0), but the second derivative vanishes along the direction 
of competition. We then have to consider the quartic term in the expansion 
of J (q) about q = q c . 

For T^O, the first and second integrals of Eq. (jSHD exist for d > (m+4)/2 
and for d > (m + 8)/2, respectively. We then have 
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(43) 



where r = (T — T c )/T c . From these asymptotic results, it is possible to obtain 
all the classical critical exponents. 

At T = 0, the first integral in Eq. (j4"Tj) exists for d > (m + 2)/2, and the 
second integral for d > (m + 6)/2. We then have 
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(44) 



where S = (g — g^)/g c - in conclusion, we have the same values for either 
classical or quantum exponents. As in the case of p ^ 1/4, the only difference 
is the critical dimension. According to an old conjecture about quantum 
critical behaviour, the quantum values of the critical exponents are given 
by the corresponding values of the classical version of the system in d + 
z dimensions, where z is a dynamical critical exponent. For m < d the 
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dynamical critical exponent is z = 1 (anisotropic case); for m = d, it assumes 
the value z = 2 (isotropic case). It should be mentioned that Hamiltonian 
formulations of the mean spherical model have been considered by some 
authors |24|. 



3.4 Decay of pair correlations 

Let us consider the system in a site-dependent field, 

H = - £ J t,t s t s t-E H t s r 

and write the partition function 
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where H (q) is the Fourier transform of H-j>, and we are omitting the vector 
notation. We then have 
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from which we obtain the pair correlations in real space, 
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As discussed by Pisani and collaborators [IB] , the analysis of (S r S r+ h), for d > 
3, below the critical temperature, leads to the introduction of a modulated 
order parameter, with characteristic oscillations for p > 1/4. A detailed 
analysis of the long-range correlations at Lifshitz point has been published 
by Frachebourg and Henkel [25]. There are also some investigations of the 
interplay between competing interactions and the decay of correlations [26J . 
We now show that these oscillations of the pair correlations in terms of 
distance are already present in the much simpler one-dimensional case. In 
the thermodynamic limit, Eq. (j4"8"|) can be written as 
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with the spherical condition 
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dq -, (50) 



27r/3 J_ w [jj, — cos q + p cos 2g] 

where p > 0, p, and /3 are written in units of J\ > 0, with the requirement 
that p > p c . Therefore, p > (1 — p) for p < 1/4, and /x > [l/(8p) + p] for 
p > 1/4. Eq. (j49p can be rewritten as 

1 / w ft+1 
(5 r , S r+h ) = — — * dw— x 2 j — , (51) 

IT flip J c W 4 - -W 3 + ^W 2 - -W + 1 



V V V 



where the contour C is the unit circle and we assuming that h > 0. The 
fourth-order polynomial in the denominator is easily factorized, with two 
roots, W\ and w 2 , inside the unit circle. After some straightforward algebra, 

we have 

, a c \ cos [fc0i + (0i - a)] .. 

(S r S r+h ) = — r exp [h In \wi\\ , (52) 

cos (0i — a) 

where a is a real phase depending on p and /i, and 

wi= |wi|exp(i0i) = — {1-A-2B), (53) 

with 



1 - 1 



A = V8p 2 -Spu+l: B= x /--2p(p-fi)- -A, (54) 

so that |wi| < 1, which leads to the well-known exponential decay. The 
spherical condition ( 1501) can be used to parametrically eliminate the chemical 
potential /j, and write the correlations in terms of T and p. In particular, 
oscillations are suppressed by 0i = 0, which is equivalent to 

/crT 1 

-k-s- P - 2 ^ (55) 

with the asymptotic value T = for p = 1/4. In Fig. 3, we draw this 
border, and indicate the region with oscillating correlations (0i ^ 0) in the 
T — p plane. From eq. fl52l) . we can write expressions for a modulation 
length, Ld = 2ii/8i, and the correlation length, £ = — 1/ln \wx\, in terms of 
temperature and the competition parameter p, which is a useful information 
to investigate the growth of modulated domains [26]. A similar behaviour 
has been found by Stephenson [27] in a calculation for the ANNNI chain. 
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Figure 3: Regions of the T — p plane of the one-dimensional model with 
simple exponential decay (6\ = 0) and oscillatory exponential decay {9\ 7^ 0) 
of the pair correlations. 



4 Conclusions 

We report an analysis of the phase diagram of a quantum mean spherical 
model in terms of temperature T, a quantum parameter g, and the ratio 
p = —J2/J1, where J\ > is a ferromagnetic interaction between first- 
neighbour sites along the d directions of a hypercubic lattice, and J 2 < 
is associated with competing antiferromagnetic interactions between second 
neighbours along m < d directions. We regain a number of results for the 
classical version of this model, including the topology of the critical line in 
the g = space, with a singular behaviour at the Lifshitz point, p = 1/4, 
for 2 < d < (m + 6)/2, which includes the case of the usual analogue of 
the Axial-Next-Nearest-neighbour Ising, or ANNNI, model. We consider in 
particular the T = phase diagram, which displays a quantum Lifshitz 
point, at p — 1/4. In the g — p phase diagram, there is a critical border, 
9c = 9c (p) for d > 2, with a singularity at the Lifshitz point if d < (m + 4)/2. 
We establish upper and lower critical dimensions and analyse the critical 
behaviour in the neighbourhood of the Lifshitz point. In one dimension, we 
derive analytic expressions for the decay of pair correlations, and determine 
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the region of modulated behaviour in the T — p phase diagrams. 
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